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Abstract 

The S-meson decay constant Jb is an important component for studying the S-meson decays, 
which can be studied through QCD sum rules. We make a detailed discussion on fs from two sum 
rules up to next-to-leading order, i.e. sum rules I and II, which are derived from the conventional 
correlator and the correlator with chiral currents respectively. It is found that these two sum 
rules are consistent with each other. The sum rules II involves less non-perturbative condensates 
as that of sum rules I, and in principle, it can be more accurate if we know the dimension-four 
gluon condensate well. It is found that /b decreases with the increment of nib, and to compare 
with the Belle experimental data on /b, both sum rules I and II prefer smaller pole 6-quark mass, 
= 4.68 lb 0.07 GeV. By varying all the input parameters within their reasonable regions and 
by adding all the uncertainties in quadrature, we obtain /b = 172]']25 MeV for sum rules I and 
Jb = 2Utll MeV for sum rules II. 
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The study of B physics plays a fundamental role both in accurately testing the Standard 
Model and in search of new physics. For such purpose, it is very import to determine 
the 5-meson transition matrix elements with high accuracy. Especially, as has already 
been emphasized in many works in the literature, the simplest matrix element J'b, which is 
defined as (0 1 (77^756! 5 (p)) = ifsPfi oi' ij^b + f^q){^\'iil5b\B{jp)) = m\fB is fundamental for 
the physics of the heavy-light quark systems, where q denotes the light u or d quark under the 
isospin symmetry. The decay constants usually are studied through their leptonic decays, 
earlier discussions of which can be found in Ref.[l|. As for Jb, its first direct measurement 
is done by Belle experiment, which shows that fs = 22913^ (stat)l'l37(syst) MeV from the 
measurement of the decay tUj. ^\. At the present, the lattice QCD calculation and the 

QCD sum rules are two suitable approaches for theoretically extracting these elements and 



their results are usual 



investigated in Refs 



y c omplementary to each other. Some typical extractions of Jb were 
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17| from lattice calculation. 



ll| from QCD sum rules and in Refs. 
The QCD sum rules approach was developed by Shifman, Vainshtein and Zakharov (SVZ) 
in the seventies of last century [isj] , and now, it becomes a widely adopted tool for studying 
hadron phenomenology. It was designed to determine the properties of the ground-state 
hadrons at low momentum transfer to the region of large momentum transfer. More explic- 
itly, hadrons are represented by certain interpolating quark currents taken at large virtuality, 
and then the correlation function (correlator) of these currents is introduced and treated in 
the framework of operator product expansion (OPE) such that the short- and the long- 
distance quark-gluon interactions are separated. The short- distance interaction is calcu- 
lated using perturbative QCD, while the long-distance one can be parameterized in terms 
of the universal non-perturbative vacuum condensates. Next, the result of QCD calcula- 
tion is matched, via dispersion relation, to a sum over hadronic states. Finally, the Borel 
transformation is introduced to suppress or even cut off the irrelevant or unknown excited 
and continuum states so as to improve the accuracy of the obtained sum rules, since little 
is known about the spectral function of the excited and continuum states. In such a way, 
the so called SVZ sum rules combines the concepts of OPE, the dispersive representation 
of correlator and the quark-hadron duality consistently that allows the derived non-excited 
hadron states' properties to be within their reasonable theoretical uncertainties. 

It is noted that one of the basic quantity for the QCD sum rules is the correlator. How to 
"design" a proper correlator for a particular case is a tricky problem. By a suitable choice of 



the correlator, one can not only obtain the right properties of the hadrons but also simplify 
the theoretical uncertainties effectively. Usually the currents adopted in the correlators are 
taken to be those with definite quantum numbers, such as those with definite J^, where J 
is the total angular momentum and P is the parity. However such kind of correlator is not 
the only choice adopted in the literature, e.g. the correlator constructed with chiral current 

to study the process B K*'y and 



261]. For different correlators, we need 



is also adopted, which is firstly proposed by Ref.jlQl 
then developed in dealing with various processes [20|- 
to check whether the theoretical estimations are consistent with each other or not, or to 
determine which correlator can lead to more accurate or more appropriate estimation for a 
particular physical process. And it is one of the purpose of the present paper, to make a 
comparative study of Jb under two typical choices for the two-point correlator. 
The first correlator is the conventional one and is defined as 



n(g2) = W d'xe'^^{0\T{q{x)j,b{x),b{0)j,qm\0) 



(1) 



and the second one is the correlator with chiral current, 

U{q')=tJ d'xe'^''{0\T{q{x){l + ^,)b{x),m{l-lMO)m • (2) 

Following the standard procedure of SVZ sum rules, we can obtain two sum rules for fs 
from the above defined correlators. 

The sum rules up to next-to-leading order (NLO) for the first correlator can be written 
as (sum rules I) 
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where mj, stands for the pole quark mass, fijR is the renormalization scale, x = m^/s and 

Cp = 4/3. M and Sq stand for the Borel parameter and the effective continuum threshold 

respectively. Since is quite small in comparison to m^ or /i/^, so terms proportional to 

mq have implicitly neglected. The function p{x) determines the spectral density of the NLO 

correction to the perturbative part and it takes the following form 
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where the dilogarithm function Li2(a:) = — Jq y ^^^^ ~ Practically, p{x) is firstly derived 



under the MS scheme, and then it is transformed to be the present form with the help of 



the well-known one loop formula for the relation between the MS 6-quark mass and the pole 

. In some references jsl. ]27|. 



quark mass, i.e. mb{piR) 



rrib 



1 + 



-4 + 3 In ■ 



it is argued that one should use 6-quark running mass other than its pole quark mass to do 
the numerical analysis. However, we think these two choices are equivalent to each other in 
principle, since we need to use their relation either to change the bound state part (to be 
determined by the pole quark mass) with the running mass or to change the hard scattering 
part (to be calculated by the running mass) with the pole quark mass. Moreover, as a cross 
check of the present obtained formulae, it is found that the sum rules I agrees with those of 
SflB6 when taH.g *e .a.e app.o.,.a«o„ and *e .a.e type of ^^ua.. ,„a.. 
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The sum rules up to NLO for the second correlator can be written as (sum rules II): 
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where s'q also is the effective threshold. Again, the terms proportional to are neglected 
due to their smallness. It can be found that the sum rules II is simpler in form than that of 
sum rules I. 

As a comparison of sum rules I and II, it can be found that by taking proper chiral 
current in the correlator, one can reduce the theoretical uncertainties to a certain degree. 
Especially, in sum rules I, the first non-perturbative term and hence the dominant contri- 
bution of the non-perturbative condensates is the dimension-three quark condensate (qq). 
And, in sum rules II, the first non-perturbative term is the dimension-four gluon condensate 
{^GG). Therefore, with proper chiral currents in the correlator, one can naturally suppress 
the non-perturbative sources and hence improve the accuracy of the obtained sum rules. 
This observation is similar to the case of light-cone sum rules for B — )■ pseudo-scalar tran- 
sition form factors where with chiral current in the correlator, one can eliminate the most 



conventional correlator (see Ref. 



uncertain twist-3 contributions 24j-|26l| in comparison to the sum rules derived by using the 



28l | for an explicit example) 



Schematically, the two sum rules I and II can be rewritten in the following simpler form 



ml Jml 

where the spectral density p*°^{s) can be read from Eqs.( l3|5i) . It is noted that the continuum 
threshold sq should be varied with the different choice of correlators. While more practi- 
cally, the Borel parameter and the continuum threshold Sq are determined in such a 
combined way that the resulting 5-meson decay constant fs does not depend too much on 
the precise values of these parameters. In addition, the following three criteria are adopted 
for numerically determining the range of of (M^, Sq): 

• Criterion (A): The continuum contribution, that is the part of the dispersive integral 
from So to oo, should not be too large, e.g. less than 30% of the total dispersive 
integral. As will be found later that this condition is well satisfied, since the resultant 
ranges for (M^,so) shall led this ratio be less than 20% for sum rules I and be less 
than 15% for sum rules II. 

• Criterion (B): The contributions from the dimension-six condensate terms shall not 
exceed 15% for Z^. 

• Criterion (C): The derivative of the logarithm of Eq.(l6]) with respect to (— 1/Af^) gives 
the _B-meson mass m^, i.e. m\ = J^2 p^°^{s)e~'^^^^^ sds j p^°^{s)e~'^^^^^ds, and we 
require its value to be full-filled with high accuracy, i.e. ^ 0.1%. 

Here, we shall not give any further constraints on sq, and we only treat it as an effective 
scale, which is roughly around the squared mass of i?-meson first excited state. 

In numerical evaluation of the sum rules I and II, we take niB = 5.279 GeV and use the 
following values for the condensates Q : 

(gg)(lGeV) = -(0.246^°:°}«GeV)^ 
(^GG) = 0.012^[]•°?^GeV^ 

TT 

{qga ■ Gg)(lGeV) = (0.8 ± 0.2)Ge1/2(gg)(lGeV). (7) 

It is found that the relative error for {^GG) is quite large. Further more, we shall use LO 
anomalous dimensions of the condensates to evaluate them up to the renormalization scale. 




FIG. 1. Allowable range of (M^, sq) for sum rules I, where mb = 4.68 GeV and the non-perturbative 
condensates are set to be theh center values. The left diagram is for criteria (A) and (B), whose 
third axis is for the ratio of the continuum contribution and the dimension-six contribution over 
the total contributions, respectively. The right diagram is for criterion (C). 




FIG. 2. Allowable range of (M^,so) for sum rules II, where nib = 4.68 GeV and the non- 
perturbative condensates are set to be their center values. The left diagram is for criteria (A) 
and (B), whose third axis is for the ratio of the continuum contribution and the dimension-six 
contribution over the total contributions, respectively. The right diagram is for criterion (C). 



A^/i? = sJ^'b ~ "^b 28|. The 6-quark pole mass is taken as: rrib ~ 4.68lo;o7 GeV, which is the 
recent world average suggested by Particle Data Group [31]. As for the running coupling 
constant, we approximate it by a one-loop form, a<((/i/_R) = 47r/('9 Indi jj^/ Annp)), where 



3l|. 



Aqcd = 0.241 GeV that is determined by reproducing as{Mz) = 0.1176 

Criteria (A, B, C) determine a set of parameters for each value of rrif,. In the following, 

we shall firstly take nib = 4.68 GeV as an explicit example to show how to determine the 

allowable ranges for and Sq and to derive Jb within these ranges. And then we shall 

present the results for several typical m?,. 

Firstly, the three-dimensional figures (fT][2|) present the allowable range of and Sq 

for sum rules I and II. The left diagram of Figs. (|Tl[2l) is for criteria (A) and (B), whose 




FIG. 3. The allowable range of (M^, sq) for sum rules I (Left diagram) and the value of /b within 
this range (Right diagram), where rrij, = 4.68 GeV and the non-perturbative condensates are set 
to be their center values. 




FIG. 4. The allowable range of (M^, sq) for sum rules II (Left diagram) and the value of Jb within 
this range (Right diagram), where nib = 4.68 GeV and the non-perturbative condensates are set 
to be their center values. 

first and second axes are for and Sq, and third axis is for the ratio of the continuum 
contribution over the total contribution for criteria (A) and the ratio of the dimension- six 
contribution over the total contribution for criteria (B) respectively. By setting a fixed value 
for So, it is found that the continue contribution increases with the increment of and the 
contribution from the dimension-six condensate term increases with the decrement of M^, 
so criterion (A) determines the upper limit of and criterion (B) determines the lower 
limit of M^. Then, a possible Borel window can be obtained from criteria (A) and (B), 
which shall be further constrained by criterion (C). Practically, the range of (M^, sq) can be 
determined numerically. The left diagram of Fig. ([3]) shows the allowable range of (M^,so) 
for sum rules I, which is obtained by sampling 500 x 400 points within G [2, 7] GeV^ and 
So G [32, 36] GeV^. And the left diagram of Fig.(jl]) shows the allowable range of (M^, so) for 
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soiGeV^) 




fsiMeV) 


rrih = 4.85 


31.7 


4.88 


119 ±6 


rrih = 4.75 


33.3 


3.52 


149 ±8 


nib = 4.68 


35.9 


2.99 


172 ± 10 


nih = 4.61 


36.9 


2.25 


193 ± 15 



TABLE I. Sum rules I for fs, where the values of sq and are those that lead to the maximum 
value of /b for a particular nih. The errors are caused by varying (M^, sq) within the region that 
is determined by setting the non-perturbative condensates to be their center values. 





soiGeV^) 




fsiMeV) 


nib = 4.85 


33.2 


2.00 


139 ± 12 


rrih = 4.75 


35.9 


2.16 


181 ± 14 


mb = 4.68 


36.9 


2.41 


214 ± 10 


rrif, = 4.61 


38.9 


2.63 


238 ± 17 



TABLE IL Sum rules II for /b, where the values of sq and are those that lead to the maximum 
value of Jb for a particular mf,. The errors are caused by varying (M^, sq) within the region that 
is determined by setting the non-perturbative condensates to be their center values. 

sum rules II, which is obtained by sampling 400 x 500 points within G [2, 6] GeV^ and 
So G [32,37] GeV^. The right diagrams of Figs. (l3P|) present the corresponding values of fs 
within its allowable region of and Sq. For fixed M^, it can be found that Jb increases 
with the increment of sq. And for fixed sq, Jb shall be steady under the reasonable region 
of the Borel parameter M^, i.e. the uncertainty is less than 3%. By varying (M^, sq) within 
their allowable region, the uncertainties of fs is less than 6% for sum rules I and less than 
5% for sum rules II. Such small uncertainty for a fixed mb agrees with the requirement of 
the reliability of sum rules. 

Secondly, we present for several typical m^, i.e. mfe=4.61 GeV, 4.68 GeV, 4.75 GeV 
and 4.85 GeV, respectively. The allowable region of (M^, sq) for each rrib can be determined 
through the same procedure described above. The results are put in TABs.l lIlini) . where the 
listed values of Sq and are those that lead to the maximum value of Jb for a particular 
rUb- And the listed errors for are determined by setting the non-perturbative condensates 
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to be their center values. By further varying the non-perturbative condensates within the 
region of Eq.Q, an extra ±1 MeV up to ±4 MeV error should be taken into consideration. 
For examples, by setting = 4.68 GeV, an extra ±2 MeV error shall be caused by those 
non-perturbative condensates for both sum rules I and II; while by setting = 4.61 GeV, 
it changes to be ±1 MeV. Furthermore, it can be found that the value of fs decreases with 
the increment of mb, and to compare with the experimental data on fs, both sum rules 
I and II prefer small 6-quark pole mass, i.e. mi, = 4.68 ± 0.07 GeV. It is found that the 
main uncertainty comes from rUb, so a better understanding of nib shall further improve the 



present sum rules. For example, in Ref. 



29| . a smaller region of rrib G [4.55,4.60] GeV is 



adopted, which can be obtained by translating the running 6-quark mass used there to the 
present pole mass, and then their result for fs is 214^5 MeV. It should be noted that such 
a small error is also due to a smaller Borel window and a smaller region of sq derived there 



29| . However if we treat (M^, sq) as correlated parameters then their separate region shall 



be broadened as shown by the left diagrams of Figs. ( 131141) . 

As a summary, we have presented a detailed discussion on fs from two sum rules, i.e. 
sum rules I and II, which are derived from the conventional correlator and the correlator 
with chiral currents respectively. Under proper parameter values, both sum rules I and II 
can lead to reasonable Jb that is consistent with the Belle experiment. At a fixed rub, the 
estimated uncertainties of fs from both sum rules I and II are less than 10%. If further 
fixing So, fs shall be quite steady versus the Borel parameter M^. Hence the contributions 
from the uncertainty sources, e.g. the continuum states and higher dimensional condensates, 
are well under control. With the help of TABs.( II|IT|) . by varying nib within a small region 
nib = 4.68 ± 0.07 GeV, we obtain = 172l| MeV for sum rules I and fe = 214^^^ MeV 
for sum rules II, where the errors from all the mentioned parameters such as nib, {M"^, Sq) 
and the non-perturbative condensates are added together in quadrature. Furthermore, the 
comparison of the sum rules I and II shows that we can improve the QCD sum rules in 
principle by using chiral current in the correlator. At the present, due to the large uncertainty 
of the dimension-four gluon condensates, such improvement on the sum rules with chiral 
current is not so clear as that of QCD light-cone sum rules for S — > P transition form 



factors 



24j-|26|. However with a more accurate values for these condensates, sum rules II 



shall become more accurate. 
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